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R. E. Hartle 
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Greenbelt, Maryland 


ABSTRACT 

The model ion-exosphere of Eviatar, Lenchek and Singer for a nonrotating 
planet with a static, centered-dipole magnetic field is generalized by permitting 
the density and temperature to vary over the baropause. In particular, the den- 
sity and/or temperature are allowed to differ at the conjugate magnetic field 
points of the baropause. For this case, the ion and electron velocity-distribution 
functions, satisfying the Vlasov equations in the drift approximation, are con- 
structed. Then, the species densities are determined from the distributions and 
compared with the corresponding results of Eviatar et. al. for the case of uni- 
form density and temperature at the baropause. In addition, the particle current 
densities, pressures and temperatures are derived. 
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I. Introduction 


A model ion-exosphere for a nonrotating planet with a static, centered- 
dipole magnetic field was constructed by Eviatar, Lenchek and Singer . 1 This 
theory has been a useful step towards a better understanding of planetary atmos- 
pheres (see, e.g., the review articles by S. J. Bauer 2 and F. L. Scarf 3 for dis- 
cussions of applications and further references to work on this subject). The 
salient features and assumptions made in the model of Eviatar et. al. are as 
follows: The base or baropause is an imaginary surface where the particle 
mean free path is of the same order of magnitude as its path length in the exo- 
sphere. The ion-exosphere is the region exterior to the baropause where 
collisions are so infrequent that they can be ignored. Within the baropause lies 
the barosphere where collisions are so frequent that the ion and electron 
velocity distributions are Maxwellian and have equal temperatures. All charged 
particles that populate the exosphere have emerged from the barosphere. Such 
complicating factors as the presence of neutral particles, ion creation and loss 
and interaction with interplanetary plasma are neglected. The following analysis 
is based on the same simplifying assumptions. 

In addition to these assumptions, Eviatar et. al. simplified the problem by 
requiring the baropause to be a spherical surface on which the density and tem- 
perature are constant. A more general baropause, of course, would consist of 
a nonspherical surface over which the density and temperature can vary. In the 
following treatment the density and temperature are allowed to vary over the 
class of baropause surfaces that are symmetric about the magnetic equator. 
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For this case, the corresponding density, pressure, temperature, etc. are de- 
rived for both species of an exosphere composed of fully ionized hydrogen. 

When the density and temperature are nonuniform on the baropause, cur- 
rents may flow from one hemisphere of the exosphere to the other leading to the 
possibility of charge buildup in the barosphere. To avoid the complication aris- 
ing from the formation of corresponding plasma sheaths at the baropause, it is 
assumed that either the currents entering or leaving the exosphere are con- 
nected in the barosphere to form closed loops or that the barosphere acts as a 
combination source and sink for current. Furthermore, only those cases of 
weak electric current are considered so that the resulting induced magnetic 
field can be ignored relative to the main dipole field. 

Following Eviatar et. al. , the particle trajectories in the ion-exosphere are 
determined by the usual nonrelativistic drift approximation. This restricts the 
exospheric models considered to those with strong magnetic fields such that the 
average Larmor radius is much smaller than the typical scale length. This also 
requires the nonmagnetic forces to be weak so that the change in particle velocity 
over a cyclotron period is small relative to the average velocity. Furthermore, 
the particle drift across magnetic field lines is assumed negligible compared 
with the average particle path length in the exosphere. 

Subject to the given boundary conditions, explicit solutions of the Vlasov 
equations in the drift approximation for the proton and electron velocity- 
distribution functions f + and f (superscripts + and - refer to protons and elec- 
trons throughout), respectively, are constructed in Section II. The electrostatic 
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field is assumed to be the usual polarization field which is then shown to be self- 
consistent with the resulting Poisson equation in Section III, where the species 
densities are derived. The corresponding longitudinal particle current densities 
and longitudinal and transverse pressures for each species are also obtained in 
Section IV. 


The cgs electromagnetic system of units is used throughout this work. 

II. Velocity- Distribution Functions: Drift Approximation 
In order to represent the boundary conditions, first consider the magnetic 
dipole field line directed from 1 to m in Figure 1. The magnetic field B is 
centered about the origin O of the planet and has an equatorial plane passing 
through the line O-m. The baropause surface, assumed continuous and sym- 
metric about the equatorial plane for simplicity, crosses the field line at points 
1 and 2 where the radial distances IrJ = |r 2 | and the latitude angles X t = - X 2 
Unless otherwise noted, subscripts 1 and 2 hereafter refer to quantities evalu- 
ated at the conjugate points ? and r 2 , respectively. The length s measured 
from point 1 in the direction of the magnetic field is given by 

^3 1/2 sin\j +Y(Xj)' 


s = (r j/^cos^Xj) 


Y(Xj) s inXj- Y(X) sinX + 3~ 1/2 ln 


^3 1/2 s in \ + Y(\) /J 


(i) 


in terms of Y (X) = (1 + 3 sin 2 X) 1/2 and the latitude angle X entering the relation 

r = r x cos 2 X/cos 2 Xj , (2) 

where r is the radial distance to the magnetic field line. Then, at r, and r 2 
the emerging distributions have the Maxwellian forms 


f 1 = / 3 * exp [- m* v 2 / 2kTj] for v • B (r x ) > 0 , 


(3a) 
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and 


(3b) 


f * = /Sj* exp [- m* v 2 /2kT 2 ] for v • B (F 2 ) < 0, 
respectively, where k is Boltzmann’s constant, m* are the component masses 
and v is the magnitude of their respective velocities v. The quantities /3.* for 
j = 1,2 equal N. (m ± /277kT.) 3/2 in terms of the baropause densities N. and tem- 
peratures T. (the superscripts + and - have been dropped since N. + = N.“ and 


The kinetic equation in the drift approximation (or guiding center approxi- 
mation) can be obtained by first introducing the drift variables (R, v^ , v M ) into 
the Vlasov equations for f*, where R is the position of the guiding center and v^ 
and V| | are the components of velocity transverse and longitudinal to the mag- 
netic field 


B = 3 



(4) 


in terms of the magnetic moment x . Then, expanding the resulting kinetic equa- 
tions in powers of the ratios rjVl, the Larmor radii to the scale length, one 
obtains the kinetic equations for f* of f* =f 0 ±+ (R b VL) f ^ + (R B ± /L) 2 f 2 ± + . . . 
which are independent of the Larmor phase. Upon neglecting terms correspond- 
ing to drift across magnetic field lines and dropping the subscript zero, the 
steady state kinetic equations for f * are 4 - 5 

3f* 1 / By* ^ m±v l 2 BB^Bf* | v ll v i BB Bf* (5) 

Bs 2B ds/Bv|| 2B Bs Bv^ 

where B is the magnitude of the magnetic field and the potential energies V* for 
the two species are given by 

V* = m* cp s ±ecp e (6) 


» 


> 
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in which e is the electronic charge and 

9 g = - GM/r (7) 

is the gravitational potential in terms of the gravitational constant G and the 

radial position r from the center of the planet of mass M. The electrostatic 
potential is taken to be the usual polarization potential of Pannekoek 6 and 
Rosseland, 7 namely <p = -(m + -m ) qp g /2e, which is shown to be self-consistent 
with Poisson’s equation in Section III. Accordingly, 

V + = V~ = (m + + rrf ) cp g /2 ■ V (8) 

and the nonmagnetic force F = - W is the same on an electron and proton. In 
this case, the drift terms ± (F x B/eB 2 ) ’Vf* ignored in Eqs. (5) can be approxi- 
mated by |F|f ± /eBL I = m + MGf ± /2r 2 eBL^ in terms of the transverse scale length L^. 
Then, the neglect of the drift terms relative to those given, say v n 3f ± /Bs~v ± f ± /L„ 
in terms of the thermal velocities v^* = (SkT* / m i ) 1/ 2 and the longitudinal scale 
L (| , is justified when the planetary models considered satisfy Bev th y|F| » 1 
or equivalently, B »1.1 x 10" l6 M/r 2 (T + ) 1/2 , assuming /L (| is of order unity. 

By the method of characteristics, 8 the general solutions of Eqs. (5) consist 
of arbitrary functions of the integration constants E* and p, 1 of the correspond- 
ing characteristic equations; that is, f* = F 1 (E*, ( u ± ), where 

E* = m 1 v 2 /2 + V, (9a) 

^ - m ± v_ l 2 /2B, (9b) 

and F* are arbitrary functions of their arguments. One notes that E 1 are the 
particle energies with v 2 = v^ 2 + Vj ( 2 and /j* are their magnetic moments or 
first adiabatic invariants. 
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The potential energies and the magnetic field are symmetric about the 

equator and their magnitudes monotonically decrease from the baropause to a 

minimum at point m. Accordingly, a particle emerging from the barosphere at 

point 1 and also passing point m will enter the barosphere at point 2 and vice 

versa. Since the particles populating the exosphere have emerged from either 

points 1 or 2, they must have E* and p* satisfying either 

a* = E* - - Vj > 0 ( 10a ) 

or 

a 2 ± ■ E* - - V 2 > 0 , (10b) 

respectively, and those particles passing point m must also satisfy 

of = E* - /^B - V > 0 , (11) 

where the quantities with subscript m are evaluated at point m. Due to the sym- 
metry of the system a* = for a given E* and ^ ; however, the separate iden- 
tification between cij* and a* is retained for generality. 

A principle noted by Aamodt and Case 9 in connection with the neutral exo- 
sphere problem is applied to construct the component distribution functions. That 
is, f* need not be continuous or single-valued functions of their arguments. With 
this in mind, consider the parts Fj* of the distribution functions f* due only to 
particles that have emerged from point 1. To insure that the orbits of these 
particles intersect point 1, F* must vanish when a* < 0. Particles passing 
point m from point 1 will not return to point 1; therefore, in segment I, between 
points 1 and m, F* must vanish when v (( < 0 and a* > 0. Accordingly, the dis- 
tribution functions Fj* in segment I satisfying these conditions and Eq. (3a) are 
F^ = fi* exp [-(E 1 -V^/kTj] S(a*) [1 - S(af) S( e 2 v n )] , < 12 ) 
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where the unit step function 


f 1 x > 0 

S (x) = •< (13) 

lO x < 0 

and e . = (-l) i+l for j = 1,2. The particles emerging from point 1 and passing 
into segment n, between points m and 2, satisfy the conditions a * > o, a * > 0 
and V|j >0. Due to the symmetry of the electric and magnetic fields, this class 
of particles with v-B>o has the same velocity distribution at conjugate mag- 
netic field points in segments I and II. Then, in segment II 


F x * = ft* exp [~(E ±- V 1 )/kT 1 ] S(a 1 ± )S(a m ± ) S( e, v„ ) . 


(14) 


Again, due to the symmetry of the system, those parts f/ of f* corresponding 
to particles emerging from point 2 are obtained for segments I and II by simply 
permuting the subscripts from 1 to 2 and 2 to 1 in Eqs. (14) and (12), respectively. 
Altogether, the total distribution functions in segment I are 

f* = yS^exp [~(E ± ~V 1 )/kT 1 ] S(a*) [l-SCo^S^v,, )] 

+ P 2 exp [-(E*- V 2 )/kT 2 ] S(a 2 ± )S(a±)S(e 2 v |1 ) , (15) 

satisfying the above requirements and the boundary condition of Eq. (3a). The 
total distribution functions in segment II are obtained by permuting the subscripts 
from 1 to 2 and 2 to 1 in Eqs. (15). At their common point m, the total distribu- 
tions for segments I and II become 


f* = i/3* exp(-Q ± /kT 1 )S( V|| ) +£* exp(-Q ± /kT 2 )S(-v || )}S(a 1 ± ) , (16) 

where Q ± = m ± v 2 /2-(V 1 -V M ) . In the symmetric case when N x = N 2 and T x = T 2 , 
the total distribution functions have the same form in segments I and II, given by 


C ± OX i /T'i it \ /i 1 n/ +\. 

I - trxp L“ - Vj)/ ki x j ) 


n\ 

v xi ; 


from which the zeroth order velocity moment is shown below to be identical to 
the density distribution derived by Eviatar et. al. 
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The question arises whether or not the kinetic Eqs. (5) are satisfied by f* of 
Eqs. (15) and the corresponding f* of segment II since they are not strictly func- 
tions of the constants E* and /a* (due to their dependance onS^ v ( | ) for j = 1,2). 
It is worth noting that a similar problem arose in connection with a solution of the 
collisionless Boltzmann equation obtained by Aamodt and Case 9 for the neutral 
exosphere. Consider the left-hand sides of Eqs. (5) which reduce to 

ifil exp [~(E ± -V 1 )/kT 1 ] exp [-(E ± -V a )/kT 2 >S(a 1 ± )S(a B ± ) 8(v„ ) ( 18 ) 

upon substitution of Eqs. (15), where S represents the Dirac delta function. As 
an aid in discussing whether or not this expression vanishes, consider the plots 
of a* and a * versus p* for v (| = 0 shown in Figure 2. The intersection point 
of a* = a* (for v,, = 0) is /a* = - (V x - V m )/(B 1 - B m ) at which point 

a* = V {1- [(B-B m )Vj/V+ (Bj-B^/Vl/tBj- Bj } . ( 19 ) 

It can be shown that the term within the curly brackets of Eqs. (19) is greater 
than zero for s x < s < s 2 along all field lines corresponding to the range 
0<kj <77/2. (To show this, one employs the expression for B obtained 
from Eq. (4) and V of Eq. (8) in conjunction with the field line coordinate 
constrained by Eq. (2)). Then, the intersection point = a m ± <0 and, referring 
to Figure 2, either S (a*) = 0 or S (a*)® 0 for /a* > - (V x - V m )/(Bj - B m ) or 
0 <iA*<-(V, - V )/(B, - B ) , respectively. Therefore, expression (18) vanishes 
and the kinetic equation is satisfied. In an analogous fashion one can show that 
f 1 for segment II also satisfies the kinetic equation. 

In general, the density and temperature vary on the baropause surface. In 
this case, the parameters Nj, N,, Tj and T 2 appearing in Eqs. (15) will be func- 
tions of position 7 which reduce to their corresponding boundary values on the 
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baropause. For example, suppose the variation of density N x is prescribed by 
some function, say N(rj) , where the tip of r x defines the baropause surface. 
Then, the function, say F(r) , replacing N x in Eqs. (15) for segment I and N x in 
the corresponding expression for segment II must satisfy F(r x ) = N(r x ) on the 
baropause. When the parameters n ,N 2 , T x and T 2 are replaced by functions 
of position, the resulting distribution functions f* will not, in general, satisfy 
the kinetic equations. However, since the kinetic equations are of zero order in 
R B VL, such a replacement is permissible when the local scale lengths of the 
functions replacing N x , N 2 , T x and T 2 are much larger than the corresponding 
Larmor radius. In this instance, the kinetic equations remain satisfied to zero 
order in R B ± /L. 

The effects of rotation, when the ion-exosphere co-rotates with the planet, 
have been treated by Melrose 10 for the case of uniform density and temperature 
at the baropause. hi this instance, the potential energies V 1 of Eqs. (6) are 
modified by the addition of the centrifugal potential energies — m* (fix r ) 2 /2| in 
terms of the angular velocity Q of the planet. It is worth noting at this point 
that the above treatment may be generalized to include such rotation when D and 
k are parallel or antiparallel. In this case, V* is symmetric about the equator 
when cp e is symmetric. However, V* will not, in general, be monotonic from 
the baropause to point m. To apply the above procedure, segments I and II must 
be divided in two parts corresponding to the points where the centrifugal force 
becomes dominant; i.e., where BVV^s vanishes. On the other hand, if cp e is not 
symmetric then V* will not be, as might be the case when the baropause radii 
| rj ^ 1 7 | (not treated here). Further segmentation may then be necessary 
since the maximum or minimum points of V* occurring at the equator in the 
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symmetric case may be shifted. In each of the above instances, the electrostatic 
potential may be determined by application of the quasi-neutrality approximation 
n + = n or by direct solution of Poisson's equation. 


III. Density Distribution 

The densities n* of the component guiding centers to zeroth order in R B */L 
are 


n* = | f 0 ± d 3 v (20) 

and, of course, are good approximations to the actual particle density when 
R b ± /L<< 1. These quadratures on f* of Eqs. (15) for segment I and the corre- 
sponding f* for segment n lead to n + = n", in which case Poisson's equation 
V 2 <p e = - 477C 2 e(n + - n* ’) is identically satisfied to the order considered. Then, 
omitting the superscripts + and -, the density ratio n/N x in segment I is 

n/Nj = exp (b n ) [1-(1 -t]) 1/2 exp (b 21 )] 


2Z ( " 1)i+1 §j eXp(b ij } {$ ( a ij>- (1 - 7 7> 1/2 exp(b 2j )$(a 4j ) 
j=l J 1 


+ O- 1) 1/2 exp (b 3 .) [0(a 2 .) -«//(a 3 .)]> 
in terms of the complementary error function 


(21) 


$ (x) = 

and a form of Dawson's integral 11 


■ ‘-(j)T 


exp (- t 2 )dt 


( 22 ) 


i/<(x)=^— J J exp ( t 2 )dt . (23) 

For j =1,2, the density ratio S. = N./N x and the arguments appearing in Eq. (21) 
are given by 
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a. j = t(V m - V)/kT .] 1/2 , a 2j =a lj/(7 _ 1) . n, (24a, b) 

a 3 j = <a?i +(V I -V„)/[(B 1 /B.-l)kT j ]}> 22 , (24c) 

a 4j = (a 2 , + (r-l)(V 1 -V m )/[(B 1 /B <i -l)kT i ]) 1 / 2 (l->))- 1/2 , (24d) 

b, , = (V, -V)AT j , b 2J = ijb./d - 17 ) , b Sj =-r»3j <24e,f,g) 

in terms of the dipole field ratios 

7} = B/Bj = (cos A.j/cos)i) 6 [(1 + 3 s in 2 X.)/(l + 3 s in 2 ^Vj)] 1/2 , (25a) 

y = B/B m = (1 + 3 s in 2 A.) 1/2 /cos 6 k. (25b) 


The density ratio n/N 1 in segment n is obtained by permuting the subscripts 1 to 
2 and 2 to 1 on the right-hand side of Eq. (21) (only the right most subscript on 
the double subscripted variables) but not on the denominator of S . = N./N x . In 
the symmetric case, when N x = N 2 and T x = T 2 , the summations in the expres- 
sions for n/Nj vanish. Thus, the density ratios n/N x in segments I and II have 
the same form given by 

n/Nj = exp (b lx ) [1 — ( 1 - exp (b 21 )] (26) 

in consonance with the result of Eviatar et. al. for this case. This relation can, 
of course, be obtained by direct integration of Eqs. (17) for f*. As noted by 
Eviatar et. al ., the term in square brackets reduces the density below that of 
the simple barometric law given by n/N x = exp(b n ). 

To exhibit some of the main differences in the density distribution along a 
field line between the symmetric model of Eq. (26) and the asymmetric model of 
Eq. (21) for segment I along with the expression corresponding to segment II, 
consider the graphs of n/N x versus latitude angle k shown in Figures 3 and 4. 

In correspondence with some of the numerical examples presented by Eviatar 
et. al. , the baropause radii r x = r 2 were chosen to be 1.1 earth radii at latitude 
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angles = \ 2 = 45° and M was taken to be the mass of earth. The solid lines 
in Figures 3 and 4 correspond to the symmetric case N x = N 2 and Tj = T 2 = 1500°K. 
To simplify the task of interpretation, only the effects of varying the density ratio 
N 2 /N x » keeping T x = T 2 = 1500°K, are considered in Figure 3. That is, the dashed 
lines in Figure 3 correspond to N 2 = O.SNj and N 2 = 2N X for Tj = T 2 = 1500 °K. On 
the other hand, in Figure 4, only the temperature T 2 is changed with respect to 
the symmetric example, where the cases T 2 = 0.5TJ = 750°K and T 2 = 2T X = 3000°K 
for Nj = N 2 are shown by dashed lines. The ratios N 2 /Nj andT 2 /Tj were chosen 
only to simplify the presentation of the results and may not be useful to any model 
of an actual exosphere. 

Relative to the symmetric case , one can observe in Figure 3 that the density 
is reduced in both segments I and H when the density is reduced at point 2 from 
the value of the symmetric example. This can be viewed as a result of the reduc- 
tion in the particle injection rate into the system. The density reduction in seg- 
ment I corresponds to a reduction in those particles that emerge from point 2 
and pass into segment I. However, at all conjugate field points, there is a greater 
reduction of density in segment II than segment I. This reflects the fact that a 
greater fraction of the particles emerging from point 2 remain in segment II than 
pass into segment I. Converse comparisons can be made when the density is in- 
creased at point 2 above the value of the symmetric example. Further, when 
N 2 = 0.5NJ , it can be noted that at point 1 (point 2) the density (n 2 ) is less 
(greater) than the density Nj (N 2 ) just below the baropause. This decrease 
(increase) at point 1 (point 2) occurs because more (fewer) particles are emitted 
at point 1 (point 2) which pass into segment II (segment I) — not returning — than 
are emitted at point 2 (point 1) which pass into segment I (segment II). A similar 
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observation can be made for the case N 2 = 2Nj . As noted by Herring and Kyle 12 
in connection with the neutral exosphere problem, such discontinuities would be 
removed in a more realistic theory where the sharp boundaries would be replaced 
by diffuse zones which account for a gradual decrease in collision frequency. 

When the temperature of the symmetric model is reduced at point 2, reduc- 
ing the particle injection rate, one observes, in Figure 4, that the total number 
of particles in the system is, of course, reduced. The decrease in temperature 
at point 2 leads to a decrease in the high energy tail of the velocity distribution 
of emergent particles with a corresponding increase in the number of low energy 
particles. This results in a relative increase in the number of particles popu- 
lating low altitudes of segment II along with a corresponding decrease at high 
altitudes. The reduction of the number of particles in the high energy tail also 
results in fewer particles being able to pass over the peak of the potential bar- 
rier at point m, leading to the density decrease in segment I. Similar compari- 
sons can be made for the case N 2 = 2Nj . As in the asymmetric examples of 
Figure 3, one observes analogous density discontinuities due, again, to the net 
flow of particles from one segment to the other. 

The above examples are not intended to represent any actual planetary ion- 
exosphere but only to indicate some of the differences between symmetric and 
asymmetric models in the simplest manner. Only one field line was considered 
because the salient differences between the two models are essentially the same 
on other field lines. 
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IV. Higher Order Velocity Moments 

The mean velocity, pressure, temperature, etc. of each component to zeroth 
order in R b Vl are functions of one or more of the velocity moments 

n<v cr p > ± = jv CT »f 0 ±d3 v ( 27 ) 

of order p= 0, 1, 2, . . . , where a corresponds to either the transverse 1 or 
longitudinal ii symbol. Referring to the distribution functions f* of Eqs. (15) 
for segment I, one notes that the velocity moments can be expressed as 

n<v CT p y = Ij 1 (p, a) - K^Cp, cr, e 2 ) + K 2 4 (p,a, e 2 ) (28) 

in terms of the quadratures 

I * (p,cr) = yS^Jv/exp [-(E ± -V 1 )/kT j ]S(a 1 ± )d 3 v (29a) 

Kj* (p.^) = P? J V/ exp [-(E* - V^/kT.] SCa^SCa^S^v,, )d 3 v (29b) 

for j = 1, 2 and -6=1,2. Permutation of all subscripts on the right-hand side of 
Eqs. (28) from 1 to 2 and 2 to 1 leads to the velocity moments in segment II. The 
first terms on the right-hand sides of Eqs. (28) can be identified as the velocity 
moments resulting in the symmetric model; i.e. , when Nj = N 2 and Tj = T 2 the 
velocity moments become Ij 1 (p.o-) since Kj (p ,a-,e^) = K 2 (p,a,e^) (the sub- 
scripts 1 and 2 are superfluous in this case). With this in mind, one may identify 
the second and third terms as those affecting the change in the moments of the 
symmetric model when N 2 4 Nj and/or T 2 / T,. An analogous identification of 
the terms of the velocity moments of segment n can be made. 

The transverse particle current densities n ^ ^ vanish to zeroth order in 
R B ± /L; however, in general, the longitudinal current densities n <y, are finite 
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to this order. The longitudinal current densities are the same at conjugate field 

points in segments I and II respectively. Accordingly, 

2 

n ^fllX = # i N i ( kT / 27m± ) 1/2 ex P ( b 4 j ) [1-(1-B m /Bj) exp(b 5 .)] ( 30 ) 

j=i 

in segments I and n, where 

b 4 i = (V,-V„)/kT jl b Sj = ^/(Bj/B^-I) (31a,b) 

are constants on a given field line. Since the cross-sectional area of a flux tube 
is proportional to l/B and constant B/ (m*) 1/2 , one notes that the par- 

ticle currents passing various points on the same flux tube are' constant; i.e., 
(d/Bs)(n<V||^>VB ) = 0. The longitudinal current densities vanish in the sym- 
metrical model and also when the density and temperature at points 1 and 2 
satisfy the constraint 


Nj/N 2 = (Tj/T 2 ) 1/2 exp (b 41 -b 42 ) 


l_( 1 _ B n ) / B 1 ) exp (b 51 ) 


(32) 


_1 - (1-B^/Bj) exp (b 52 )_ 

It is worth noting that this expression is simply Nj/N 2 = (T x /T 2 ) 1/2 to zeroth 
order in (V x -V m )/kT. when | (V x -V m )/kT. | << 1 for j = 1, 2. 


When there is a net particle flow from one hemisphere to the other, there 
will always be an electric current since<V || )>~/ < (y 1| X = (m + /nf ) 1/2 . These currents 
must be weak in the sense that the induced field B' satisfies I B'| / |b| < < 1. From 
Ampere's law, the magnitude of the induced field B' can be approximated in 
terms of its scale length L' by B'~47ren| <y ( j ^>“ | L ' . Then, requiring B'/B<<1 
and using the propertyn <y M y /B = constant, one has 

L'/Bj « l/(477en| <v, (33) 
which relates a necessary restriction on the ion-exospheric models considered. 
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For the asymmetric examples considered in the previous section, the par- 
ticle current densities at are: 2.6 x lo 6 !^ and -3.4 x lo 6 !*! cm -2 sec -1 for 
T 2 = 0.5^ = 750°K and T 2 = 2T X = 3000°K when = N 2 ; 2.0 x 10% and 
-3.9xl0 6 N 1 cm‘ 2 sec‘ 1 for N 2 = 0.5^ and N 2 = 2N X when Tj = T 2 = 1500°K. 
When a conservative value n <V M ^> - ~ 5 x 10 6 N X is assumed L'/Bj N x << 10 12 cm 4 
gauss -1 . Then, for these cases, the class of models that can be treated by the 
above are restricted by BjNj >> 10 -3 gauss cm -3 when a conservative value 
L' ^ 10 9 cm (the approximate length of a field line) is assumed. 


To zeroth order in R VL , the transverse and longitudinal pressures Pj* and 


p,p respectively, are defined in terms of the velocity moments by 

p l = j m± f v i { o ± d 3 v = Y m±n < v i y 
Pi* =m± |( v ir<V||> ± ) 2f 0 ±d3v =m±n [ < v n 2 >* - (< v ,i>V] 

and, for the symmetric case, with N = N 2 and T x = T 2 , they become 


Pi = NjkTj exp (b n ) 


1_(1_77> 1/2 (l + -J 1 ?-b 21 ^ exp (b 21 ) 


Pil = NjkTj exp (b u ) [1_(1 -t 7> 3/2 exp (b 21 )] , 


(34a) 

(34b) 


(35a) 

(35b) 


explicitly, in both segments I and n. The superscripts have been omitted since 
the species pressures are equal on the symmetric model. The transverse and 
longitudinal pressures for the more general case Nj / N 2 and/or Tj / T 2 are 
given in the Appendix for completeness. Since the difference (Pn -Pj_ ) on the 
symmetric model can be expressed as 


Pil “Pi = N i kT i (l- 7 ]) 172 



exp (b 21 / V) . 


(36) 
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one notes that p, , > at all points s x < s < s 2 and thatpu = p^ = NjkTj at points 
1 and 2, their isotropic value just below the baropause. Furthermore, one can 
show that both the transverse and longitudinal pressures decrease, along a flux 
tube, from their values at the baropause to a minimum at the equator. 


The transverse and longitudinal temperatures and T M , respectively, fixed 
by p^ = nkT^ and p| ( = nkT ( , , can be written as 


Tj. =Tj 


1 - 




1/2 


(i^- b 2i) 


exp (b 21 ) 


1-(1-7]) 1/2 exp(b 21 ) 


(37a) 


T ll = T i 


1 _( 1 _t 7 ) 3/2 exp (b 21 ) 
1_(1_T)) 1/2 exp (b 21 )_ 


(37b) 


on combining Eqs. (26) and (35). It is clear that is less than its value T x at 
the baropause. On the other hand, T (| increases along a flux tube from Tj at the 
baropause to a maximum at the equator . 
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Appendix 


The transverse and longitudinal pressures on the assymetric model are ob- 
tained by carrying out the operations of Eqs. (34a) and (34b), respectively, on f + 
of Eqs. (15) for segment I and the corresponding f ± for segment n. Accordingly, 
for segment I, one has 


Pl = NjkTj exp (bjj) 


1 “ ( 1 - 7 ?) 1/2 U+TT^-b-J exp (b 21 ) 


1 


2lJ F v 21 ; 


\ J2 (-lV + 1 N.kT. ^(a,.) 

(* + 7^- b 2j) $ ( a 4j) +7T “ 1/27 7 a 4j ex P(- a 4j 2 ) 


- ( 1 -'l) 1/2 exp (b, ) 


+ (7- 1) 1/2 exp (b, ) 


, 1 + T y " b 3i) ^ (a 2j> ~ ^ < a 3 j >} 


- tt~ 1/2 y {a 2j exp (a 2 2 ) - a 3j exp (a 3 2 . )} 


3 j 


(Al) 


Pll = NjkTj exp ( b n ) [ 1 -( 1 - tj ) 3/2 exp ( b 21 )] - m *!! (( v ,,),*) 2 

+ I Z + (^) 1/2 a U ex P(- a i?> 


- (l-7]) 3/2 exp (b ) 


- (7- 1) 3/2 exp(b 3j ) 


/4 \l/2 

$ ( a 4j) + (-J a 4j exp(-a 4 2 ) 


/ 4 \ 1/2 

’/< a 2 j) -V'( a 3j ) -( — j a 2 j eX P< a 2 j 2 ) 


/ 4 \ 1/2 

+ M a 3j exp(a 3J ) 
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where the superscripts have been omitted since the species pressures are equal 
(one notes, on combining Eqs. (21) and (30), that the second term of Eq. (A2) 
is independent of particle type). The corresponding pressures in segment II are 
obtained by permuting the subscripts 1 to 2 and 2 to 1 on the right-hand sides of 
Eqs. (Al) and (A2) (only the right most subscript on the double subscripted 
variables). 
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Figure 2-The arguments and aj 1 of the step functions versus jjr when V|| 
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Figure 3-Density ratio n/Nj versus latitude angle k along a magnetic field line crossing 
the baropause at 45° latitude and 1. 1 earth radii. The solid line corresponds to Nj = 
N 2 . For comparison, the dashed lines correspond to N 2 = 0.5N 1 and N 2 = 2N r The 
baropause temperatures T, - T, = 1500 °K in each case. 
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Figure 4-Density ratio n/Nj versus latitude angle k along a magnetic field line crossing 
the baropause at 45° latitude and 1.1 earth radii. The solid line corresponds to T l ~ 
T 2 13 1 500°K. For comparison, the dashed lines correspond toT 2 ”0.5T 1 ~ 750 Kand 
T 2 - 2Tj - 3000°K. The baropause densities N 1 ~ N 2 in each case. 




